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Abstract. We derive the exact expression of the diffusion coefficient of a self-gravitating Brownian gas 
in two dimensions. Our formula generalizes the usual Einstein relation for a free Brownian motion to the 
context of two-dimensional gravity. We show the existence of a critical temperature T c at which the diffusion 
coefficient vanishes. For T < T c , the diffusion coefficient is negative and the gas undergoes gravitational 
collapse. This leads to the formation of a Dirac peak concentrating the whole mass in a finite time. We 
also stress that the critical temperature T c is different from the collapse temperature T» at which the 
partition function diverges. These quantities differ by a factor 1 — l/N where N is the number of particles 
in the system. We provide clear evidence of this difference by explicitly solving the case N = 2. We also 
mention the analogy with the chemotactic aggregation of bacteria in biology, the formation of "atoms" in a 
two-dimensional (2D) plasma and the formation of dipoles or "supervortices" in 2D point vortex dynamics. 

PACS. 05.45.-a Nonlinear dynamics and nonlinear dynamical systems 



1 Introduction 



Recently, there was a renewed interest for systems with 
long-range interactions pQ. These systems have a very 
strange thermodynamics and exhibit peculiar features that 
are very different from those of more familiar systems with 
short-range interactions like neutral gases and plasmas. 
One striking property of systems with unshielded attrac- 
tive long-range interactions is their ability to self-organize 
spontaneously into large-scale coherent structures. Some 
examples are provided by stars, globular clusters and galax- 
ies in astrophysics, jets and vortices (e.g. Jupiter's great 
red spot, gulf stream,...) in two-dimensional geophysical 
flows, bacterial aggregates in biology, clusters in the Hamil- 
tonian Mean Field (HMF) model etc. As a result of spatial 
inhomogeneity and non-extensivity, the ordinary thermo- 
dynamic limit (N — > +oo with N/V fixed) is clearly irrel- 
evant for these systems and must be reconsidered [2J3J. On 
the other hand, the statistical ensembles are generically in- 
equivalent and the choice of the relevant ensemble must be 
addressed specifically. Therefore, systems with long-range 
interactions are special and escape the ordinary rules of 
thermodynamics. Among all types of systems with long- 
range interactions, self-gravitating systems are probably 
the most fundamental 4J. These systems have a complex 
thermodynamics and present interesting phase transitions 
between "gaseous" and "clustered" states [5] . They can in- 
deed undergo catastrophic collapse below a critical energy 
E c in the microcanonical ensemble [6] or below a critical 



temperature T c in the canonical ensemble [7] when gravi- 
tational attraction overcomes diffusive effects. 

In astrophysics, the dynamics of self-gravitating sys- 
tems is basically described by the Newton equations where 
the acceleration of a particle is equal to the gravitational 
force (per unit of mass) created by the other particles 
8 9 10 . These equations have a Hamiltonian structure. 
For a large number of particles N ^> 1, wc cannot fol- 
low the motion of each particle in detail and we must 
have recourse to statistical mechanics |4I5] , An isolated 
self-gravitating system is conservative (evolving at fixed 
energy E) so that the proper statistical ensemble is the 
microcanonical ensemble. This is the correct description 
of stellar systems such as globular clusters and galaxies. 
The Hamiltonian A-stars model has a very long history 
starting with Newton's Principia Mathematica in 1687. 

In a recent series of papers, Chavanis & Sire [11-21] 
have introduced and systematically studied a model of 
sclf-gravitating Brownian particles. In this model, the par- 
ticles interact gravitationally but they also experience a 
friction force and a stochastic force which mimick a cou- 
pling with a thermal bath of non gravitational origin. 
Therefore, the basic equations of motion consist in a set 
of A coupled Langevin equations. This system is dissipa- 
tive (evolving at fixed temperature T) so that the proper 
statistical ensemble to consider is the canonical ensemble. 
This Brownian model could find applications in the pro- 
cess of planetesimal formation in the solar nebula [22]. In 
this context, the dust particles experience a friction with 
the stochastic force due to turbulence and, when 
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the dust layer becomes dense enough, self-gravity must 
be taken into account. On the other hand, some interest- 
ing analogies have been found between a self-gravitating 
Brownian gas and the process of chemotaxis in biology 
|17l23j , the formation of large-scale vortices in 2D turbu- 
lence |24l25j and the Bose-Einstein condensation in the 
canonical ensemble |26j . At a more academic level, Brow- 
nian motion is a fundamental process in physics (pio- 
neered by the works of Einstein and Smoluchowski) and 
it is clearly of interest to investigate the situation where 
N Brownian particles are coupled by a long-range po- 
tential of interaction, like the gravitational interaction. 
These various arguments may justify the study of the self- 
gravitating Brownian gas model. 

Random walkers in interaction described by N coupled 
Langevin equations are also studied in soft matter physics 
to compute transport properties of systems consisting in 
many interacting particles such as supercooled liquids or 
the dynamics of colloids in solution [27128) . In these ex- 
amples, the potential of interaction is short-range and the 
system is spatially homogeneous at equilibrium. An inter- 
esting problem consists in determining the effective diffu- 
sion coefficient of a particle of the system. This is a very 
difficult problem and essentially perturbative approaches 
have been developed. For example, Dean & Lefevre 29J 
consider a regime of weak coupling (small (3 with fixed 
(3p) where diagrammatic expansions can be carried out. 
Interestingly, their approach suggests that the system will 
experience a glassy behaviour (identified by the divergence 
of the relaxation time and the vanishing of the diffusion co- 
efficient) at some critical temperature T c . However, since 
the expansions are valid at high temperatures, it is not 
clear whether their extrapolation to lower temperatures 
of order T c is justified. 

In this paper we show that, for the self-gravitating 
Brownian gas in two dimensions, the diffusion coefficient 
of a particle can be calculated exactly for any temperature. 
It is given by the following expression 

where T c is the critical temperature 

Cm 2 

k B T c = (N-l) — . (2) 

For T 3> T c , when self-gravity becomes negligible, we re- 
cover the Einstein relation [30]. In that case, the parti- 
cles have a diffusive motion (corresponding to an evapo- 
ration of the system) slightly modified by self-gravity. For 
T = T c , the diffusion coefficient vanishes and for T < T c 
it becomes negative implying finite time collapse. In that 
case, the system forms a Dirac peak containing the whole 
mass in a finite time. These different regimes have been 
studied in detail in the mean field approximation by ana- 
lytically solving the Smoluchowski-Poisson system |12ll9j . 
Here, we present complementary results that are valid be- 
yond the mean field approximation. We stress that, con- 
trary to the case of systems with short-range interactions 



like those evocated above [27I28|29] . we calculate the dif- 
fusion coefficient of a particle in a gravitational system 
that is out-of-equilibrium and spatially inhomogeneous. 
The context is therefore very different from the case of 
soft matter physics [27|28I29 ], We also stress that the re- 
sult (HI)-© only holds in two-dimensional gravity. Grav- 
ity is known to be critical in two dimensions because, di- 
mensionally, the gravitational potential u = Gra 2 ln(rjj ) 
does not depend on the distance (the distance enters in 
a dimensionless logarithmic factor). This is the intrinsic 
reason why the Virial of the gravitational force and the 
diffusion coefficient can be calculated exactly in 2D grav- 
ity. 

The paper is organized as follows. In Sec.[2l we develop 
a many-body theory of Brownian particles in interaction 
and derive exact kinetic equations valid for an arbitrary 
binary potential of interaction (Sec. 12. ip . In Sec. 12. 3i we 
derive the general expression of the Virial theorem for a 
Brownian gas with two-body interactions. When we con- 
sider the gravitational potential in two dimensions (Sec. 
12. 4|) . the Virial theorem takes a very simple form from 
which we can deduce the exact expression of the diffusion 
coefficient (Sec. l2.5| ). The mean field approximation is con- 
sidered in Sec. 12.61 In Sec.[3j we obtain the exact equation 
of state of a self-gravitating gas in two dimensions. This 
equation of state, first derived in plasma physics for elec- 
tric charges [31 , is well-known. It is usually derived in 
the canonical ensemble from the partition function |32l4j . 
Here, we show that the same expression can be obtained 
in the microcanonical ensemble from the density of states. 
We also generalize its expression by allowing the parti- 
cles to have different masses (or charges). In Sec. HJ we 
investigate the existence of statistical equilibrium for a 
two-dimensional self-gravitating system in the canonical 
ensemble. We show that the partition function exists only 
above a temperature T* . This temperature differs from the 
critical temperature T c appearing in the Virial theorem 
and in the equation of state by a factor 1 — 1/N where 
N is the number of particles in the system. We provide 
clear evidence of this difference by explicitly solving the 
case N = 2 in Sec. [5] This difference has been overlooked 
in the literature. However, this is essentially a curiosity 
because, in the thermodynamic limit N — ► +oo, the two 
temperatures coincide. 

2 Many-body theory 
2.1 Exact kinetic equations 

In a space of dimension d, we consider a system of N 
Brownian particles confined within a domain (box) of vol- 
ume V and interacting via a binary potential |2|3j . Their 
dynamics is described by the coupled Langevin equations 



= -£«v Q V a U(n,...,r N ) + y/2D^R a (t), 

(4) 
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where £ Q is a friction coefficient, D' a a diffusion coeffi- 
cient and R Q (i) a white noise satisfying (R a (t)) = and 
(Ra t i(t)Rpj(i?)) = 8 itj 5 a> p5(t - t') where a = 1, ...,N la- 
bels the particles and i = 1, d the coordinates of space. 
The particles interact via the potential 

C/(n, ...,rjv) = ^2 m am u(\r a - r p \), (5) 

a<f3 

where u(\r — r'|) is a binary potential of interaction de- 
pending only on the absolute distance between the parti- 
cles. The TV-body Fokker-Planck equation describing the 
evolution of this system is 



dP, 



N 



N 



Of 



£ v. 



a=l 



N fl 

Y — 



D 



OPn 
dr a 

, dPN 



N 



£ a PNV a 



(6) 



where Pjv(ri, vi, r/v, vjv, t) is the iV-body distribution 
function and F Q = — ^— V a U is the force by unit of mass 
experienced by particle a. In order to obtain the canonical 
distribution at statistical equilibrium, we must impose the 
Einstein relation 



D' 



£, a k B T 



(7) 



This shows that the temperature T is a measure of the 
strength of the stochastic force in Eq. ([4}. Then, the sta- 
tionary solution of the Fokker-Planck equation ([5]), can- 
celling independently [3] the advection term (l.h.s.) and 
the "collision" term (r.h.s.), is the canonical distribution 



Pv(ri, Vi, ...,r N , vjv) 



-/3H(ri,vi,...,rjv,vjv) 



, (8) 



where H is the Hamiltonian 

N v 2 

H = 2_j m a~^~ + m a mpu(\v a — vp\). (9) 

a—l a<p 

In the overdamped limit where £ Q — > +oo, we can ne- 
glect the inertia of the particles in Eq. ((4]) and we get 



dt 



-ti a V a U(v 1 ,...,r N ) + ^/2DoR. a {t), (10) 



where we have introduced the mobility \i a = l/(£ Q m Q ) 
and the spatial diffusion coefficient D a = D' a /(^. The N- 
body Fokker-Planck equation describing the evolution of 
this system is 



dP N 
dt 



N a 

y— 

^ dr a 



dP, 



N 



dr a 



Ha"NT{ 



(11) 



where P/v(ri, ...,r^,t) is the iV-body distribution in con- 
figuration space. In order to obtain the canonical distri- 



bution at statistical equilibrium, we must impose the Ein- 
stein relation Q 



fi a 1 



(3. 



(12) 



Using the expression of the mobility, the Einstein relation 
can be rewritten 



D„ 



k B T 



(13) 



This expression can also be obtained from Eq. ([7]). Then, 
the stationary solution of the Fokker-Planck equation (fTTj) 
is the configurational part of the canonical distribution 



P/v(ri, ...,rjv) 



f 



Z(J3) 



-/3(7(r 1 ,...,r JV ) 



(14) 



If we introduce the one and two-body probability distri- 
butions 



P a {r,t) = J P N {r x , ...,r N ) dr 7 , 



(15) 



P a ,p(r,r',t) = J P N (r u ...,r N ) J] dr 7 , (16) 

we find from Eq. (fTTj) that the one-body distribution func- 
tion ([15]) satisfies a kinetic equation of the form 



dP a _ d_ 
dt dr 



x J dr'P ai/3 (r,r',t)^(|r-r'|) 



(17) 



An alternative derivation of this equation is given in Ap- 
pendix El This equation is exact and takes into account 
statistical correlations encapsulated in the two-body dis- 
tribution (|16p . As a result, this equation is not closed since 
it involves a distribution of higher order. The complete hi- 
erarchy of equations for the reduced distributions is given 
in 12151. 



2.2 The potential energy tensor 

We introduce the potential energy tensor 

Wij = — ^2 m a m P J drdr' 



Oil 



xP Q>/3 (r,r',^ — (|r-r'|). 



(18) 



For a multi-components system, a necessary condition for 
the Fokker-Planck equation to admit a stationary solution 
is that the ratio fi a /D a be independent on a. This ratio is 
identified with the inverse temperature /3. 

2 This is the usual iV-body potential energy tensor [8] aver- 
aged over the noise. 
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Since u(\r — r'|) depends only on the absolute distance 
between the particles we get 



Wij — — m a mp I dvdr 1 



u'(|r-r'|) 



xP a ^(r,r',t)xi— —{xj-x'A. (19) 

|r — r'| 

Interchanging the dummy variables a, (3 and r, r', we have 
equivalently 



Wij = m a m0 / drdv' 



xP^r'^x f'^ r/|) (^-^.). (20) 
|r — r'| 

Summing Eqs. (fT9| and (f20|) . we obtain 
u'(\r - r'\) 

x^(r,r',t) " _ U (Xi - 4)(xj - x'j). (21) 

Under this form, the potential energy tensor is manifestly 
symmetric 

Wij = Wji. (22) 
The trace of the potential energy tensor is 

Wu = - ^ ^ m a mp J drdv' 



xP a ^(r,r',ty(|r-r'|)|r-r'| 



(23) 



2.3 The exact Virial theorem 

We introduce the tensor of inertia 

hjit) = J Pa(r,t)£ a m a XiXjdr, 

ru J 



(24) 



where £ = jj J2 a is the average friction coefficient. 
From Eq. (|T7|) . we obtain 



dP a 
dx k 



£,iij = - ^ / dr(xiS k j +XjS k i) 

dr'P a ^(r,r',t)^— (|r - r'|) 



This expression incorporates the friction coefficients £ a . 
This generalization proves to be necessary in order to obtain a 
closed expression of the Virial theorem for a multi-components 
system. For identical particles, the expression (|24[) coincides 
with the usual tensor of inertia of a iV-body system [8] aver- 
aged over the noise. 



where we have used an integration by parts (the boundary 
term cancels out since the current of diffusion vanishes 
on the box due to the conservation of the normalization 
condition J P a dr = 1). On the other hand, 



dP a 

{x l 5 k j + Xj6ki)-Q^-dr 



2Sij — <j> P a (xidSj + XjdSi), 



(26) 



and 



~ X] / dr ( x i § kj + XjSki) ^ m a i 

/Ou 
dr'P a , l3 {r,r',t) — {\r-r'\) = 2W ij . (27) 

We thus obtain 

U,. = 2Nk B TS lJ + 211,, 
- k sT I P a {x % dS + XjdSi), (28) 

which is the general expression of the Virial theorem for a 
Brownian gas of particles in interaction in the overdamped 
limit. Introducing the local pressure 



(29) 



the Virial theorem can be rewritten 

£4 = 2Nk B T6 l3 + 2Wij 

p(xidSj + XjdSi). (30) 

The scalar Virial theorem is obtained by contracting the 
indices leading to 



1 ^ 



where 



£7 = dNk B T + W u - (t pr ■ dS 



I(t) = P a (r,t)£ a m a r 2 dr, 



(31) 



(32) 



is the generalized moment of inertia (including friction 
coefficients) . If the pressure is constant on the boundary 
of the domain, we have 

<j> pr ■ dS — pi, <j> r • dS = pb J V • r dr = dpbV. (33) 



(25) More generally, we introduce the notation 



p =dvV pr ' dS > 



(34) 



which can be identified with a kinetic pressure. Then, the 
scalar Virial theorem can be written 



1 



= dNk B T + W vl - dPV. 



(35) 
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At equilibrium (I — 0), we have 

dNk B T + Wu - dPV = 0. 



(36) 



In the absence of interaction, we recover the perfect gas 
law 



PV = Nk R T. 



(37) 



2.4 The gravitational potential 

The gravitational potential in d dimensions is given by 



<0 



1 G 

' d-2^ d ~ 2 



u(0 = Gln£ (d = 2), 



(38) 



(39) 



where G is the gravity constant (whose value depends on 
the dimension of space). The gravitational force is given 

by 



AO 



G 



(40) 



Inserting this expression in Eq. (|2 1[) , we find that the grav- 
itational potential energy tensor reads 

Q ^ p 

Wij = -— ^ m a m,f3 I drdr' 
1 ^ J 

n I i ,\ ( X i ~ X 'i)( X j ~ X 'j) 
XP «./3( r ' r .*) | r _ r ,|d 

and that the trace of the potential energy tensor is 

Q p 

Wu = -— ^ m a m/3 / drdr' 
1 «*fi J 



(41) 



xP a ,/?(ry,t) _ a . 



For d ^ 2, we have 



W ti = (d- 2)W, 



(42) 



(43) 



where W is the usual potential energy [8] averaged over 
the noise. In that case, the scalar Virial theorem (|3"5|) can 
be written 



1 



£J = dNk B T + (d-2)W- dPV. 



(44) 



For d = 2 a "miracle" occurs and we get the simple exact 
result 



G 



E 



m a mp, 



(45) 



where we recall that J2a^/3 = Yl a =i S/3^ Q - We note that 
this expression, contrary to Eq. (|43j) , does not depend on 



the configuration of the system but only on the masses of 
the particles. We define a critical temperature 



kaT r = 



47V 



For equal mass particles, we have 



k B T c = (N-1) 



Gm 2 



(46) 



(47) 



Then, the scalar Virial theorem ([33]) in d = 2 can be 
written 



yi = Nk B (T-T c )-PV. 



(48) 



At equilibrium, in a bounded domain, we get the exact 
equation of state 



PV = Nk B (T-T c 



(49) 



where P is the kinetic pressure (fM]) . In Sec. [31 we shall 
obtain this equation of state directly from the partition 
function and deduce, by identification, that P also rep- 
resents the thermodynamical pressure. Since P > 0, we 
conclude that a necessary condition for the system to be 
at statistical equilibrium is that T > T c . For T < T c , 
we cannot have statistical equilibrium since the pressure 
would be negative. In fact, for T < T c , we have I < e < 
so that I{t) — at a finite time t = t en d implying that 
the system forms a Dirac peak containing all the mass 
in a finite time. Therefore, statistical equilibrium is not 
possible for T < T c . However, the condition T > T c does 
not guarantee statistical equilbrium. We shall see in Sec. 
|4]that statistical equilibrium is possible only for T > 
where is strictly larger than T c . For T c < T < T„, the 
system collapses (or has a peculiar temporal behaviour) so 
that Eq. (|49]) . corresponding to I = 0, is never satisfied. 

In an infinite domain, assuming that initially the par- 
ticles arc in a confined region of space, the pressure at 
infinity vanishes (P — 0) and the scalar Virial theorem 
reduces to 



(50) 



-^I = Nk B {T-T c ). 

This equation is readily integrated leading to 
m = ^Nk B{T _ Tc)t + m _ 



(51) 



For T < T c , the moment of inertia vanishes (I = 0) at 



Lend 



ANk B (T c — T) ' 



(52) 



implying that the system forms a Dirac peak containing 
the whole mass in a finite time t enc i- For T > T c , I(t) — * 
+oo for t — > +oo indicating that the system evaporates. 
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2.5 The diffusion coefficient in d = 2 



2.6 The mean field approximation 



In this section, we restrict ourselves to a one component 
system. In that case, 



pr dr, 



(53) 



is the usual moment of inertia with p — NmP a for any 
a = 1, N. The mean squared displacement of a particle 
(any) is given by 



2 = Jpr 2 dr = _I_ 

x ' M M 

Thus, from Eq. ([50]) . we obtain 

d(r 2 ) _ Ak 



dt 



-r^(T-T c ), 



(54) 



(55) 



where T c is given by Eq. (|4"T)) . After integration, we get 
4fc B 



(r 2 ) = -^(T-T c )t + (r 2 ) . 



(56) 



The diffusion coefficient of the particle is defined by 

(r 2 ) ~ 4£>(T)i, (t-^+oo). (57) 



Therefore, it is given by the exact expression 



D{T) 



(T-T c ). 



(58) 



In the absence of gravity (G = T c = 0), or for high temper- 
atures (T T c ) where gravitational attraction becomes 
negligible with respect to thermal motion, we recover the 
usual expression of the diffusion coefficient given by the 
Einstein relation 



D = 



(59) 



However, for smaller temperatures, gravitational effects 
come into play and the expression of the diffusion coeffi- 
cient is modified. Interestingly, there exists a critical tem- 
perature T c at which the diffusion coefficient vanishes El- 
For T > T c the diffusion coefficient is positive so that the 
system evaporates. For T < T c the diffusion coefficient 
becomes negative implying finite time collapse to a Dirac 
peak containing the whole mass in a time 



4fc s (T c -T)' 



This time behaves like (T c — T)^ 1 for any T < T c and 
diverges at the critical point T c . At T = T c the Dirac 
peak is formed in infinite time. These different regimes 
have been studied in [12116119] by solving the (mean field) 
Smoluchowski-Poisson system in two dimensions. 

4 This is physically different from the vanishing of the diffu- 
sion coefficient in the case of colloids 29 which is due to the 
close packing of the particles (steric hindrance) in the glassy 
phase (the particles cannot move freely) while in the present 
situation we have a collapse due to an attractive potential. 



The preceding results are exact for self-gravitating Brow- 
nian particles in d = 2 dimensions, whatever the number 
N of constituents. In particular, they take into account 
statistical correlations. However, in most works on self- 
gravitating systems, one usually considers a mean field 
approximation where the equations of the problem can 
be simplified. In general, this approximation is valid for 
N ^> 1. In this section, we briefly describe how the mean 
field approximation can be implemented and how the re- 
sults are modified. 

The essence of the mean field approximation is to as- 
sume that the two-body probability distribution can be 
written as the product of two one-body probability distri- 
butions according to 



P a , / 3(r,r',t)=P a {r,t)P f3 (r',t). 



(61) 



This approximation allows one to close the hierarchy of 
equations at the level of Eq. (TP?)) . For systems with long- 
range interactions, it can be shown that this approxima- 
tion is exact in a proper thermodynamic limit N — » +oo 
with i] = fiNu* and e = E/(u«N 2 ) fixed, where m* is the 
typical value of the binary potential. For the gravitational 
potential in d dimensions where = Gm 2 /R d ~ 2 , the 
thermodynamic limit corresponds to N — > +oo in such a 
way that the dimensionless temperature r\ = (3GMm/ R d ~ 2 
and the dimensionless energy e = ER d ~ 2 /GM 2 are fixed. 
One can always rescale the quantities of the problem so 
that the coupling constant (G in gravity) scales like it* ~ 
1/N while E - N, T - 1 and V ~ 1 rj. In this limit, the 
factorization (fBTjl is valid up to terms of order 1/N. Since 
we consider a large number limit, we can also extend the 
sum in Eq. pT|) over all the particles. This yields 



dP a _ d 
dt dr 



D a 



dP a 
dr 



p a m a P a (r,t) 



'^Y. m P J dr'P p (r',t)u{\r-r'\) 



If we introduce the density 

p(r,t) = ^P«(r,<)r 



(62) 



(63) 



(60) W e get 



dP a _ d_ 
dt dr 



dP a 

D a —^- + p a m a P a (r,t) 
dr 



xV J dr'p(r',t)u(\r - r'| 



(64) 



This equation can be rewritten in the form of a mean field 
Smoluchowski equation 



dP a 
dt 



d_ 
dr 



dP 

D a — — + p, a m a P a {r,t)V<P 
dr 



(65) 
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where the potential is determined by the density according 
to 



#(r,i)= / p(r',t)u(\r - r'\)dr' 



(66) 



For the gravitational potential, the preceding equation is 
equivalent to the Poisson equation 



A<P = S d G P . 



(67) 



Therefore, in the mean field approximation, we have to 
solve the multi-components Smoluchowski-Poisson system 
l(55 p -(|S7 p . The steady state of Eq. (JBSJ) is the mean field 
Boltzmann distribution 



P„(r) = A Q e- /3m °* (r) . 



(68) 



The single component Smoluchowski-Poisson system has 
been studied in [11-17,19-21] and the two components 
Smoluchowski-Poisson system has been studied in [TB] . 

To establish the expression of the Virial theorem in the 
mean field approximation, we can use a procedure similar 
to the one developed in Sec. 12.31 The only change is the 
factorization (|6ip and the replacement of X)/3^ Q by ^2p- 
Thus, the previous relations remain valid provided that 
Wij is replaced by 



where $(r,t) is given by Eq. (|6"B")) . In particular, for the 
gravitational interaction we have 



^ X i X 'j| X j X „■ ^ 

xp(r,t)p(r ,t) — — J —, 



and the trace of the potential energy tensor is 



W 



MF 



For d ^ 2, 



TT", 



G 



MF 



I 



drd? 



,p(r,i)p(r',i) 



r — r 



i\d-2 



(d-2)W M *, 



where W Mt is the mean field potential energy 



W MF = - j ptdr, 



1 



with 



#(r,t) 



G 



d-2J r-r 1 



rdr'. 



For d = 2, we get 
w mf = _G 



J drdr'p(r,t)p(r',t) = - 



GM 



f2 



(70) 



(71) 



(72) 



(73) 



(74) 



(75) 



Therefore, in the mean field approximation, the scalar 
Virial theorem of a self-gravitating Brownian gas in d = 2 
can be written as in Eq. (|48|) with the critical temperature 



k R T MF = 



GM 2 
AN 



For equal mass particles it reduces to 

MF GNm 2 
k B T™ F = ■ — . 



(76) 



(77) 



These expressions can be directly obtained from Eqs. (|46|) 
and dHJ) by replacing J2 a ^/3 h Y T, a in E q- & or 
by replacing TV — 1 by N in Eq. (|47|) since the mean field 
approximation is valid for N ^> 1. We note that the mean 
field results are relatively close to the exact results even 
for a moderate number of particles. Equation (|48|) is al- 
ways valid and finite N effects just slightly shift the critical 
temperature T c . Since N— 1 = N(l — 1/iV), the correction 
is of order 1/N, which is precisely the domain of validity 
of the factorization hypothesis (161|) as shown in [2|3j at 
a more general level. This corroborates the observation 
that the mean field approximation provides a good de- 
scription of systems with long-range interactions such as 
self-gravitating systems. 



(69) 3 The exact equation of state 



In this section, we derive the exact equation of state of 
a self-gravitating gas in two dimensions. We extend the 
derivation given by Salzberg |32j and Padmanabhan [4] 
in two respects: 1. we consider a multi-components sys- 
tem while the previous authors assume that the particles 
have the same mass 2. we treat both the canonical and 
the microcanonical ensembles while the previous authors 
only consider the canonical ensemble. For comparison, we 
also discuss the case of a two-dimensional plasma made of 
electric charges [3T] . 



3.1 Canonical approach 

For the gravitational interaction in two dimensions 



U(rx, ...,r N ) = G mjrnj In|r.j - r 

i<j 



(78) 



the configurational part of the partition function in the 
canonical ensemble is given by 



Z{fi,V)= / e 



; 111 \Ti—Tj | 



\{dv k . (79) 



fc=i 

To avoid the divergence of the partition function at large 
distances, we assume that the system is enclosed within 
a box of radius R. The following calculations also assume 
implicitly that the partition function converges at small 
distances. The existence of statistical equilibrium states 
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will be discussed in Sec. [4] Using a trick due to Salzberg With the transformation y = r/R, it can be written 
|32j . we set y = r/R. Then, the partition function can be 

rewritten ' g(E, R) = R 2N g{E', 1), 



Z(P,V) = R 2N e~ 0GlnR ^ 



with 



^ fe=l 



(80) 



E' = E- GlnR^2', 



l<3 



where the last integral is now independent on R. From this The entropy is defined by 
expression, we find that 



dZ _ 2N 
dR~~R 



Z(/3,R). 



The thermodynamic pressure is defined by 
1 dlnZ 



P 



P dV ' 



V = ttR 2 . 



S{E,R) = k B \ng(E,R). 
(81) According to Eq. (|89|) , we have 

S(E,R) = 2Nk B lixR + S(E',l). 
The temperature and the pressure are given by 



(82) 



From Eqs. (|82[) and (fSTj) . we obtain the exact equation of 
state of a two-dimensional multi-components self-gravitating 
gas 



1 
T 



dS 
dE 



P = T 



N,V 



dVj 



(89) 
(90) 

(91) 
(92) 

(93) 



N,E 



p = ^L 

(3V 



(3G^ 



i<3 



with V = ttR 2 . From Eqs. ([92]), ((90j) and ([93])-a we find 
that 



This can be written 

PV = Nk B {T-T c ), 
with the critical temperature 



x>y m l m J 



4N 



(83) 



(84) 



(85) 



dS_ 

dR 



N,E 



2Nk L 
R 



G ^ — > 

— — > TOiTO,-. 

TR ^— ' 



(94) 



Using Eq. (|M)) -b. we obtain 



PV = Nk B 



T{E) 



G 



2Nk, 



■E 



irnmj 



(95) 



For equal mass particles, we have 



k B T c = (N-l) 



Gm 2 



This returns the expression (|84[) obtained in the canoni- 
cal ensemble with the critical temperature (|55|) and 
Therefore, the equations of state coincide in the two en- 
(86) sembles for any number of particles. This result was not 
^ obvious at first sights. For two-dimensional self-gravitating 

This returns the equation of state (|49|) obtained by the systems, the microcanonical and canonical ensembles are 
kinetic approach. Note that for a single species system in equivalent at the thermodynamic limit with N — > +oo 



the mean field approximation, the equation of state 



but they are not equivalent for finite values of N. For ex- 



can also be obtained by solving the Boltzmann-Poisson ample, we shall explicitly show in Sec. that the caloric 
equation and computing the pressure at the edge of the curves for N — 2 calculated in the microcanonical and 



box 



3.2 Microcanonical approach 

The Hamiltonian of a self-gravitating system in two di 
mensions is 



canonical ensembles differ. Yet, the equation of state is 
the same for any N. 



3.3 The case of electric charges 



The case of electric charges is obtained by taking G = — 1 
" i and by making the substitution m <-> q. Therefore, the 

H = 7: m i v i + G 2_. rnifrij In |rj — rj-|. (87) energy of interaction in two dimensions reads 



N 



i<j 



The density of states in the microcanonical ensemble is 
defined by 



N 

g{E,R) = / 8(E - H) JJ dr k dv k 

"* k=l 



U(ri,...,r N ) = -^ftgj-lnlri-r,-! 



(96) 



Kj 



5 This is because the series of equilibria of a two-dimensional 
self-gravitating gas does not present turning points contrary to 
3D self-gravitating systems; see, e.g., [12] , 
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The previous results remain valid for electric charges in- 
stead of point masses, provided that we make the above 
mentioned substitutions. In particular, Eqs. |83|) and (|95j) 
are now replaced by 



of N particles). Finally, it is easy to determine the criti- 
cal temperature of a non neutral 2D plasma consisting in 
iV + charges +e and iV_ charges — e. From, Eqs. (|97|) and 
(fTU0|) . we find 



P = 



Nk B T 
V 



A 

2N 



(97) 



knTc. = — 



1 - 



N 



(103) 



Using the condition of electroneutrality 



N 



4 Existence of statistical equilibrium in the 
canonical ensemble 



= 0, 



(98) 



we find that 




i<j 



2^ 



For a finite two-dimensional self-gravitating system, there 
exists statistical equilibrium states for any value of the 
energy in the microcanonical ensemble (Hamiltonian sys- 
tems) . Indeed, if the system is enclosed within a box, the 
density of states g{E) is convergent for any E. By con- 
trast, in the canonical ensemble (Brownian systems) equi- 
librium states exist only for sufficiently high temperatures 
T > . Although the existence of a collapse temperature 
T* is well-known, there is some ambiguity in the litera- 
ture concerning its precise value. In particular, we will 

_ , . . „ . .. . . show that the collapse temperature T* does not exactly 

Ihcrcforc, the equation of state oi a neutral two-dimensional ... i , , m ■ i j 

, ... coincide with the critical temperature l c introduced pre- 

plasma can be written viously. We will also show that the collapse temperature 

T* is difficult to determine for a multi-components system 
while it takes a simple expression when the particles have 



1^ 

»=i 



(99) 



PV = Nk B (T - T c 



(100) 



with the critical temperature 



the same mass 



T N a 2 
AN 



(101) 4.1 Statistical equilibrium state for T > Ti 



If the plasma consists in N/2 charges +e and N/2 charges 
— e, we get 



k B T c = — . 



(102) 



These results have been first derived by Salzberg & Prager 
[3Tj in the canonical ensemble. The calculations of Sec. 13.21 
show that they can also be derived in the microcanonical 
ensemble. Since P > 0, a necessary condition for the sys- 
tem to be at statistical equilibrium is that T >T C . How- 
ever, this condition does not guarantee statistical equilib- 
rium. We shall see in the sequel that statistical equilibrium 
exists only for T>T* = 2T C . For T < T* there is no equi- 
librium state and the system collapses. This leads to the 
formation of N/2 pairs (+, — ) corresponding to non inter- 
acting "atoms" . By contrast, for T > T* the system is fully 
ionized. For comparison, a gas of self-gravitating Brown- 
ian particles collapses to a single Dirac peak containing 
the N particles for T < = jtzjT c while it remains dif- 
fuse for T > . Since the collapse in plasma physics leads 
to the formation of individual pairs, we understand qual- 
itatively why the collapse temperature k B T* = e 2 /2 does 
not depend on N (it corresponds to the condition to form 
one pair). By contrast, since the collapse in gravity leads 
to the formation of a single Dirac peak containing all the 
mass, the collapse temperature T* = NGm 2 /4 depends 
on N (it corresponds to the condition to form a cluster 



In this section, using the arithmetic-geometric mean in- 
equality, we show that the partition function Z(JS) is con- 
vergent for T > T\. We extend the method developed by 
Kiessling 34j to the case of a multi-components system. 
First, we note that 

e ~i 3G J2 i< j ™i m 3 111 l r i- r il 

i 

-| 1/N 



n 



Now, by the arithmetic-geometric mean inequality 



l/N 
i ' 



(104) 



(105) 



we get 



j_V^ -f/3G£ jV .m i m i ln| ri -r.,-| 



- PGrriirrij In |r* — rj \ 



(106) 
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From this inequality, we find that the configurational part 
of the partition function satisfies 



^E/^n/v 



- (3Gmirrij In \ vi — rj \ 



(107) 



The integrand is maximum for = yielding 

Z< -^Y^ J dri IiJ e-^ (3Gm ^ lnx dx 

En /v^ Gm * ro ' ^ 

Ell C x x -^^ Gm ^dx. (108) 

.• jj: Jo 



itR 2 

2n 2 R 2 r R 



If 



NGm imj 



(109) 



for all i, j with i j, then all the integrals converge and 
the partition function is finite. Therefore, a sufficient con- 
dition for the existence of statistical equilibrium is that 
T > Ti with 



NGmimu 



(110) 



where I and II denote the two most massive particles in 
the system. In that case, we obtain the following bound 
on the partition function 



Z < 



(111) 



For equal mass particles, a sufficient condition for the ex- 
istence of statistical equilibrium is that T > Ti with 



k B Ti = 



(112) 



In that case, we recover the bound on the partition func- 
tion given by Kiessling [34 : 



Z < ttR 2 



1 - ^l3Gm 2 



N-l 



(113) 



4.2 Collapse for T <T 2 

We now show that the partition function is divergent for 
T < T 2 . We first rewrite the configurational partition func- 
tion in the form 



Z = J dr x ...dr N J[ 



1 



(114) 



We select a particle (allowed to move over the whole box) 
and approach the N — 1 other particles at a distance e of 
the first [34135) . The contribution of this configuration to 
the partition function behaves like 



Z ~ R 2 e^ N -V \{ - 



•i \ y hSGmimA-ZSN-l) 



(115) 



Considering now the limit e — > 0, we see that the partition 
function diverges if T < T 2 with 



k B T 2 



(116) 



(117) 



4(JV- 1) 

For equal mass particles, we find that 
k B T 2 = N^f. 

4.3 The collapse temperature T* 



In conclusion, the partition function converges for T > Ti 
and diverges for T < T 2 . For T 2 < T < T±, we have no re- 
sult in the general case. However, for equal mass particles, 
we have T\ = T 2 = T* with 



k B T« = N 



Gm 2 



(118) 



Therefore, the partition function converges for T > 
and diverges for T < T*. For T > T«, we have a diffuse 
gas and for T < T„ the most probable distribution is a 
Dirac peak containing all the particles. We note that T* 
differs by a factor 1 — 1/N from the critical temperature 
T c appearing in the equation of state ([54)) . Indeed, 



N 



N - 1 



(119) 



This point has been overlooked in the literature. In Sec. 
O we shall provide a clear evidence of this difference by 
considering explicity the case N — 2 where the partition 
function can be calculated exactly. Finally, we note that 
the pressure ([54]) at the collapse temperature T* is equal 
to 



PV = fc B TL 



(T = T»). 



(120) 



This corresponds to the ideal pressure due to TV* = 1 
effective particle (the Dirac peak) containing the whole 
mass when the system has collapsed. 

For a multi-components system, the collapse temper- 
ature T» seems difficult to obtain. The case of a two- 
components system has been considered by |18j in the 
mean field approximation and an explicit expression of 
the collapse temperature has been obtained in that case. 
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4.4 The case of 2D plasmas 



5.1 The microcanonical ensemble 



For self-gravitating systems, the partition function diverges 
for T < T* due to the collapse of all the particles at the 
same point. Thus, below T*, the most probable structure 
is a Dirac peak containing the whole mass. By using an 
argument similar to that of Eq. (| 1 1 5|) . we see that if we 
concentrate only a fraction of the particles (e.g. a pair), 
we obtain a smaller collapse temperature, so that this con- 
figuration is less favorable. 

By contrast, for neutral plasmas with N/2 charges +e 
and N/2 charges — e, the divergence of the partition func- 
tion corresponds to the formation (collapse) of individual 
pairs (+,—)■ Therefore, we can obtain the collapse tem- 
perature by considering the statistical mechanics of only 
one pair (N — 2). This is done it the next section and the 
analysis yields 



The Hamiltonian of a system of two particles with mass 
mi and m,2 in gravitational interaction can be written 



H 



1 



-MV 2 



1 



-uv 



Gm lnr, 



(124) 



where (R, V) denote the position and the velocity of the 
centre of mass and (r, v) the position and the velocity of 
the reduced particle. On the other hand 



M = mi + 7772 , 



fj, 



77717772 

mi + fni 



(125) 



denote, respectively, the total mass of the particles and 
the mass of the reduced particle. We shall first compute 
the hypersurface of phase space with energy less than E, 



knT> 



(121) 



Therefore, the partition function of the iV-body system 
converges for T > T* and diverges for T < T*. For T > 
T*, we have a fully ionized gas and for T < T» the most 
probable distribution is N/2 pairs (+,—). We note that 
T» differs by a factor 2 from the critical temperature T c 
appearing in the equation of state (|100[) . Indeed, 



(122) 



Since the equilibrium states exists only for T > T», the 
equation of state (|100[) is valid only above T* . For T < 
(including the region T c < T < ) , the partition function 
diverges due to the collapse of pairs of point particles of 
opposite sign. Finally, we note that the pressure (jlOOp at 
the collapse temperature is equal to 



PV = ^Nk B T* 



(T = T, 



(123) 



This is the pressure created by an ideal gas of 2V* = N/2 
non-interacting pairs (+, — ). Finally, for T < T c , the Virial 
theorem shows that the N* pairs collapse to a single point 
in a finite time (if the motion of the charges in the plasma 
is described by Eq. (fit)]) ') while we have no result for T c < 
T < T*. In that range of temperatures we probably have 
a homogeneous distribution of pairs. 



r(E) = [ dRdVdrdv. 

Jh<e 



(126) 



Noting x 2 = {M/2) 1 / 2 Vi for i = 1,2 and x> = (///2) 1/2 ^ 
for i = 3, 4, the preceding quantity can be rewritten 



r(E) 



AttR 2 
HM 



dr 



d A x. (127) 



| x|| <\/ E— Gm\m2 lnr 

The last integral represents the volume of a four-dimensional 
hypersphere with radius \JE — Gm\m2 lnr. Therefore, we 
obtain 



, s 8ir 2 R 2 

nm = — rpVi 



(E — Gui\ m-i lnr) rdr, 



(128) 



where r m denotes the maximum distance accessible to the 
reduced particle and V4 is the volume of a four-dimensional 
hypersphere with unit radius. From the general expression 



V d 



T d/2 



r{d/2 + \y 



(129) 



we find that V4 — it 2 /2. The density of states in the mi- 
crocanonical ensemble is given by g(E) — dr/dE. Using 
Eq. (I128|) . we obtain 



5 The case N = 2 

The statistical mechanics of two particles in gravitational 
interaction (binary star) was considered by Padmanabhan 
[1] in d = 3 dimensions. He showed that this "toy model" 
exhibits phenomena that are representative of more realis- 
tic stellar systems with a large number of particles. Here, 
we shall extend his analysis in d = 2 dimensions (the one- 
dimensional case is treated in Appendix [Cj . For N = 2 
particles, it is possible to compute the density of states 
and the partition function exactly. This will allow us to 
illustrate the preceding results on an explicit example. 



&ir 4 R 2 
fj,M 



(E — Gmi7772 In r) rdr. (130) 



The range of integration r m is such that E—Gm\m,2 In r > 
and r < R. Therefore, r m = e E/G mi m 2 if E < Q mini2 m R 
and r m = R if E > Gm\vri2 lni?. Introducing 



A = 8tt 4 G, e 



E 



G?77l 7772 



t = 



k B T 

G?77i 7772 



we obtain 



g(e) = AR 2 



In r) rdr, 



(131) 



(132) 
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Fig. 1. Caloric curve in microcanonical (MCE) and canoni- 
cal (CE) ensembles for the binary star model N = 2 in two- 
dimensions. The caloric curves are similar in the two ensembles 
(without negative specific heat region) so that there is no phase 
transition contrary to the situation in d = 3 [415] . There exists 
equilibrium states for all energies in MCE and for tempera- 
tures t > t* = 1/2 in CE. It can be of interest to compare 
these curves with the caloric curve obtained within the mean 
field approximation valid for N 3> 1 which presents a similar 
behaviour [33112] , 



with r m = e c if e < lni? and r m = R if e > lni?. The 
entropy is given by 



5(e) = k B \ng(e), 

and the temperature by 

1 _ dS_ 1 _ d{S/k B ) 
T ~ dE ~* 1 ~ d~e ' 

Finally, the pressure is given by 
,<95 



P = T 



dV 



V = irR- 2 



(133) 



(134) 



(135) 



The caloric curve t — t(e) and the equation of state 
in the microcanonical ensemble can be obtained from the 
exact expression (|132|) of the density of states. For e < 
lni?, we obtain 



7(e) = 



AR 2 



GmiTO2 



t = 



PV = k B T. 
For e > In R, we obtain 

rs ARi ( in 1 
9(e) = — U-lnR + - 



(136) 

(137) 
(138) 

(139) 



with 



t = e-lnR+ -. 

2 



PV = 2k B (T-T c ), 

, Gmim 2 
k B T c = . 



(140) 
(141) 

(142) 



The caloric curve t(e) in the microcanonical ensemble is 
represented in Fig. [1] 



5.2 The canonical ensemble 

For N = 2 particles in gravitational interaction, the par- 
tition function in the canonical ensemble is given by 



Z = 



-f3M\- e -pH^r e -f3mim 2 G In r dRdY drdv. 



Integrating on R and on the velocities, we get 



Z = 



f3 2 mirti2 J 



r l-0Gm 1 m 2 dr 



(143) 



(144) 



The partition function is finite if, and only if, T > T* with 



For equal mass particles, we have 
, T Gm 2 



(145) 



(146) 



These expressions agree with the general results of Sec. HI 
For TV = 2, we have T% — T 2 — even if the particles 
have different masses m\ ^ m,2- For T > T„, the partition 
function is explicitly given by 



Z = 



8ir 4 



J^4—0Gm 1 m2 



fi 2 m\mi 2 — (3Gmim,2 



(147) 



This expression corresponds to the bound obtained in Eq. 
(jlll|) . The pressure can be computed from Eqs. ([52")) and 
(fT47|) leading to 



with 



PV = 2k B (T-T c ), 

Gm 1 m 2 
k B T c = . 



For equal mass particles, we get 

Gm 2 



knT r 



(148) 



(149) 



(150) 



These expressions are consistent with Eqs. (p?4" |) - (]86p . We 
also check explicitly that T* 7^ T c (for N — 2, we have 
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Fig. 2. Equation of state for the binary star model N = 2 in 
two-dimensions. Equilibrium states exist only for t > t* = 1/2 
so that the portion of the curve t c — 1/4 < t < t* = 1/2 is not 
physical. Note that the pressure does not vanish at t = t* = 
1/2. It corresponds to the pressure due to the Dirac peak (one 
effective particle) that forms at this critical temperature (see 
Sec. US} , 



0.4 




Fig. 3. Distribution of energies for the binary star model N = 
2 in two-dimensions in the canonical ensemble. We have plotted 
this distribution for different values of the temperature. 



= 2T C ). Therefore, T c is smaller than the collapse tem- 
perature below which the partition function diverges. 
The curve P(T) is plotted in Fig. H 

The caloric curve (e)(t) in the canonical ensemble can 
be obtained from the exact expression (|147[) of the parti- 
tion function. The average energy is given by 



(E) 



d\nZ 

d(3 



,d\nZ 



Using Eq. (I147|) , we obtain 

{At - 3)t 



2t- 1 



■InR. 



(151) 



(152) 



The caloric curve (e) (t) in the canonical ensemble is repre- 
sented in Fig. [T] and is compared with the caloric curve in 
the microcanonical ensemble. We note that, for N = 2, the 
specific heats C = d(E)/dT diverges like oc (T — X 1 *) -2 at 
T = T*. By contrast, in the limit N — > +oo where a mean 
field approximation can be implemented, we find from Eq. 
(40) of Sire & Chavanis [12] that C oc (T - T*)- 1 . 

Finally, the distribution of energies at temperature T 
in the canonical ensemble is given by 



P(E) 



1 



Z(T) 



9{E)e 



-0E 



(153) 



Using P(E)dE — P(e)de, this can be rewritten 



P(e) 



1 



/(e)e- 



e/t 



Gmim.2 Z(t)' 
Using Eqs. (fT36|) . (fl39)l and (fT47l) . we find that 
2t- 1 



P(e) = 



4t 3 



,(2-l/t)e 



(e < 0), 



(154) 



(155) 



P(e) 



2t-l 



(2e+l) e - £/ *, (e>0). (156) 



We have taken R = 1 to simplify the expressions. The 
distribution of energies is represented in Fig.[3]for different 
temperatures. 

Taking G = — 1 and m; = ft , these results also describe 
a two-dimensional plasma with 1 charge +e and 1 charge 
— e. According to Eq. (|145p . the partition function exists 
only for T > T* with 



= — , 



(157) 



while the critical temperature appearing in the equation 
of state (|148p is given by 



kaTr 



e 

T 



Thus 



(158) 



(159) 



More general studies [36] show that the expression (|15T|) 
of the collapse temperature remains valid for N/2 charges 
+e and N/2 charges — e where N is arbitrary. Indeed, the 
"catastrophic collapse" at T = T» corresponds to the for- 
mation of N* — N/2 non-interacting pairs (+,— ) such 
as those studied here individually. It can therefore be ob- 
tained by studying the collapse of just one pair. There- 
fore, in 2D plasma physics, = 2T C for any N (for a 
neutral plasma with two components). By comparison, in 
2D gravity = jf^Tc (for a one component gas). 
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5.3 Relevance of a statistical description for N = 2 
particles 

We may wonder about the relevance of studying a system 
of N = 2 particles using statistical ensembles. 

For Brownian particles, the iV-body distribution func- 
tion is governed by the Fokker-Planck equation (O whose 
steady state is the canonical distribution J5J). This is valid 
for any N. For N = 1, we recover the ordinary Brown- 
ian motion of a free particle submitted to a friction and 
a noise. The equilibrium state is the Maxwell-Boltzmann 
distribution predicted by statistical mechanics. For N = 2, 
we have two Brownian particles interacting via self-gravity 
and this system can be studied using the canonical ensem- 
ble. More generally, for Brownian particles, the canonical 
ensemble is justified for any N. For Hamiltonian systems, 
the relevance of a statistical description for a small number 
of particles is less clear. For example, N = 2 particles in 
gravitational interaction have a simple deterministic mo- 
tion satisfying the Kepler laws. A statistical description 
may be justified, however, when the particles are hard 
spheres enclosed within a box so that their motion is very 
irregular due to collisions between themselves or against 
the walls (see the discussion in the review of Padmanab- 
han [4]). In that case, the microcanonical ensemble may 
be employed. 

In any case, the density of states and the partition 
function can be defined mathematically for any N and 
it is interesting to compute these quantities exactly for 
N = 2 and compare these results with the ones obtained 
for N — > +oo by using a mean field approximation. For 
3D self-gravitating systems, it is shown in 4 5 that the 
caloric curves (in microcanonical and canonical ensembles) 
obtained with N — 2 particles are already very close to 
those obtained for a large number of particles. Similarly, 
we have shown in the present paper that the same obser- 
vation holds in d — 2 and d = 1 dimensions. This means 
that the mean field results are reliable even for a small 
number of particles. 

6 Conclusion: analogies with chemotaxis and 
point vortex dynamics 

In this paper, we have obtained the exact expression of the 
diffusion coefficient of a gas of self-gravitating Brownian 
particles in two dimensions. Two-dimensional gravity is a 
rare example where exact results can be obtained 0. The 
diffusion coefficient exhibits a critical temperature T c be- 
low which it becomes negative indicating finite time blow- 
up. In a finite domain, this critical temperature is slightly 
smaller than the collapse temperature T* below which 
the partition function diverges. However, this difference 
is small, the ratio being of the order of 1 — 1/N, and the 
two temperatures coincide at the thermodynamic (mean 

6 Of course, the results of this paper can be easily extended 
to power-law potentials of the form = G/£ d ~ a in d di- 

mensions. In that case, the critical dimension is d cr it = a + 1. 



field) limit N — > +oo. We may wonder what happens dy- 
namically in the region T c <T < T* since there is no equi- 
librium state in that case while we cannot prove finite time 
blow up from the Virial theorem. The system collapses for 
T < but it is possible that, for T c < T < T*, it takes an 
infinite time to form the Dirac peak. We stress that the 
computation of the diffusion coefficient brings only partial 
information on the dynamics of the self-gravitating Brow- 
nian gas. For example, it does not give any information 
concerning the precise evolution of the density distribution 
of the particles. To obtain the dynamical evolution of the 
density profile p(r, £), we need to solve the iV-body equa- 
tions (JSJ-jl]) or make a mean field approximation, valid 
in the proper thermodynamic limit N — > +oo, where the 
problem reduces to the study of the Smoluchowski-Poisson 
(SP) system. Therefore, the complete description of the 
self-gravitating Brownian gas remains highly complicated 
and requires in general some approximations even if an 
exact result (H}-© can be obtained. 

The results obtained in this paper can also be rele- 
vant in mathematical biology [37]. They apply to a sim- 
plified version of the Keller-Segel model [35] modelling 
the chemotaxis of bacterial populations. The Keller-Segel 
model can be derived from a system of coupled Langevin 
equations of the form [39140] : 

= X V a c + V2DR a (t), (160) 



^ = -kc + D c Ac + h^2S(r-r a (t)). (161) 

a — 1 

Here, r a (t) represents the position of a biological entity 
(cell, amoeba, bacteria,...) and c(r,t) is the concentration 
of the secreted chemical. The cells have a diffusion motion 
with diffusion coefficient D and they also move along a 
positive gradient of chemical (attractive chemotaxis \ > 
0). The chemical is produced by the cells with a strength 
h. It is in addition degraded at a rate k and diffuses with 
a diffusion coefficient D c . The usual Keller-Segel model is 
recovered from these equations by making a mean field 
approximation leading to 

^-=DAp- X V .(pVc), (162) 



— = -kc + D c Ac + hp. (163) 

Some authors have considered a simplified chemotactic 
model where the equation (|161[) for the chemical is re- 
placed by a Poisson equation 

N 

Ac=-X^2mS(r-r a (t)). (164) 

This can be justified in a limit of high diffusivity of the 
chemical and, in addition, for sufficiently dense systems 
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|41|40j . In that case, the chemotactic problem becomes 
isomorphic to the study of self-gravitating Brownian parti- 
cles described by Eq. (jlOp . with a proper re- interpretation 
of the parameters. In particular, the concentration — c of 
the chemical plays the role of the gravitational potential 
<P. We have the additional correspondanccs A «-> SaG and 
X [im. All the results of this paper can then be extended 
to the biological context provided that we make the ap- 
propriate substitutions. Furthermore, the consideration of 
the dimension d — 2 is particularly justified in biology 
since cellular organisms are often ascribed to move on a 
plane. In chemotaxis, we rather use the mass as a control 
parameter instead of the temperature. For example, the 
critical mass corresponding to the critical temperature T c 
is 



M r 



8nD 



(165) 



and the collapse mass corresponding to the collapse tem- 
perature T» is 



8nD 



(166) 



For a large number of cells (M ^> m) the two values co- 
incide. For M < M*, a box-confined system reaches an 
equilibrium state. For M > M*, the system collapses and 
forms a Dirac peak corresponding to chemotactic aggrega- 
tion. The adaptation of our results to the biological con- 
text is straightforward [12] so it will not be discussed in 
more detail here. 

Let us finally consider the case of point vortices in two- 
dimensional hydrodynamics (see, e.g., [24]). The Hamilto- 
nian of the point vortex gas is Q 



(167) 



where the coordinates (x, y) are canonically conjugate [43] . 
The partition function in the canonical ensemble is given 

by 



Z = 



■ S l<7 T*Ti ln l r i- r al 



N 



(168) 



fe=i 



We note that there is no kinetic (quadratic) term in the 
Hamiltonian (|167|) in the usual sense, so that the temper- 
ature can take both positive or negative values [13] , By 
contrast, only positive temperatures are allowed in plasma 
physics and gravity. Let us first assume that the point vor- 
tices have the same circulation 7. At negative tempera- 
tures, the partition function is formally equivalent to that 
for self-gravitating systems. Therefore, at very negative 



7 For brevity, we do not take into account the contribution of 
the images (in a bounded domain) because most of the follow- 
ing results are obtained by considering configurations where 
the vortices form compact clusters. They are therefore inde- 
pendent on boundary effects. 



inverse temperatures, point vortices have the tendency to 
"attract" each other and group themselves in a single ag- 
gregate (supervortex) of circulation A7. This corresponds 
to a regime of high energies. This is similar to the for- 
mation of a Dirac peak in the gravitational problem for 
T < . Thus, the partition function exists if and only if 
> /3* with 



0* 



8tt 

'ivy 



(169) 



This is the exact equivalent of the collapse temperature 
(|118p as can be seen by taking G = — 1/(2%) and m 7. 
At positive temperatures, the partition function is for- 
mally equivalent to that for a non-neutral plasma of equal 
charges. Therefore, at large positive inverse temperatures, 
point vortices have the tendency to "repell" each other 
and accumulate on the boundary of the domain. We now 
consider N vortices with circulations +7 and N vortices 
with circulations —7. At positive temperatures, the parti- 
tion function is formally equivalent to that for a neutral 
plasma. Therefore, at large positive inverse temperatures, 
the vortices of opposite sign have the tendency to "at- 
tract" each other and form N dipoles (+,—). This corre- 
sponds to a regime of high negative energies. This is sim- 
ilar to the formation of N pairs (+, — ) in plasma physics 
for T < T*. Thus, the partition function exists if and only 
if < #, with 



0, = %. 

7 



(170) 



This is the exact equivalent of the collapse temperature 
(|121|) . At negative temperatures, the vortices of same sign 
have the tendency to "attract" each other and form an ag- 
gregate of positive circulation A7 and an aggregate of neg- 
ative circulation — A 7 (vortex dipole). This corresponds 
to a regime of high positive energies. This is similar to the 
formation of a Dirac peak in the gravitational problem, 
except that we have here two peaks: one with positive cir- 
culation and one with negative circulation. The partition 
function exists if and only if > /3* where 0„ is given 
by Eq. (|169[) like in the case with one cluster of N parti- 
cles. Therefore, the point vortex gas is very rich because it 
exhibits features similar to self-gravitating systems at neg- 
ative temperatures and features similar to plasma systems 
at positive temperatures. We emphasize, however, that 
point vortices form a Hamiltonian system described by the 
microcanonical ensemble so that the control parameter is 
the energy (not the temperature) and the quantity of in- 
terest is the density of states (not the partition function) . 
It can be shown that the density of states converges for 
any accessible energy. Therefore, the above collapse tem- 
peratures (|169p - (|170p represent lower and upper bounds 
on the caloric curve 0(E) which correspond to E — * +00 
and E — > — 00 respectively. The question to know what 
happens for < — -^j- or > ^ has a priori no sense in 
vortex dynamics since we cannot impose the temperature 
of the vortex gas. 
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A Alternative derivation of the kinetic 
equation 

In this Appendix, we provide an alternative derivation of 
the exact kinetic equation Eq. (|17p . The Langevin equa- 
tion (I10|) governing the time evolution of particle a in the 
strong friction limit can be written 

dr, 



=-m at x a V& a >(r a ) + ^/2D^R a {t), (171) 



where 



^)(r,t) = ^m^(|r-r^(t)|), (172) 

denotes the exact potential created by the other particles. 
We define the one and two-body probability distributions 

by 



P a (r,t) = («J(r-r a (t))>, 



(173) 



F Qi/3 (r,r',t) = (S(r - r a (t))6(r> - r p (t))). (174) 
Differentiating Eq. (I173| with respect to time, we obtain 
dP a 



dt 



-V.(r tt (i)5(r-r a (t))). 



(175) 



Substituting from Eq. (TrTTj) in Eq. (TITS]) , we find that 

Pip 

^ = -V • (v^Ra(t)i(r - r a (t))) 
at 

+V • (m af i a V& a \r a (t))6(r - r a (t))). (176) 

The first term is the standard term that appears in deriv- 
ing the Fokker-Planck equation for a pure random walk; 
it leads to a term proportional to the Laplacian of P a : 



- V ■ (V2D^n a (t)5{r ~ r a (t))) = D a AP a 

To evaluate the second term, we first write 

*< a >(r,t)= f dr'^m^dr-r'D^r'-r^W). 



(177) 



(178) 



This yields 

(V#<«>(r a (t))5(r-r a (t))) = 
/ dr' ]T m/jVuflr - r'|)(<5(r - r a (i))<5(r' - r^t))) 

= Jdr'J2 m p P a ,p(r, r', *)Vu(|r - r'|). (179) 



Substituting Eqs. l(T77|) and ([179]) in Eq. ([T76]) . we get the 
exact kinetic equation (I17| . 

We can easily extend this approach so as to take into 
account the inertia of the particles. In that case, the dy- 
namical evolution of particle a is described by stochastic 
equations of the form 



dr a 
dt 



(180) 

V$ {a Hr a (t)) + ^2D^ a {t). (181) 



We define the one-body probability distributions in phase 
space by 



P a (r,v,t) = (<f(r-r a (t))a(v-v a (t))>, 



(182) 



where the brackets denote an average over the noise. Simi- 
larly, the two-body probability distribution in phase space 
is 

P Qi(3 (r, v; r', v', t) = (S(r - r a (t))5(v - v„(i)) 

5(r>-r (t))5(-v'-Mt)))' (183) 
Taking the time derivative of P a , we get 
BP r) 

-gf = -fr- <**(*)<*(' - ra(*))5(v - v a (t))) 

• (v a (t)«(r - r a (t))£(v - v a (t))>. (184) 

av 



Inserting the equations of motion (|180|) - (|18ip in Eq. I|184|) . 
we obtain 



at 



9v 



<9r 



(v a (t) <J(r-r Q (i))<y(v-v Q (t))) 



c9v 



(£ a v a (t) tf(r-r a (t))*(v-v a (t))) 



(V^(r Q (t)) 5(r - r a (t))5(v - v a (t))) 



(V2D^R Q (t) «(r-r a (t))«(v-v (i))). (185) 



The first two terms are straightforward to evaluate. The 
fourth term is the standard term that appears in deriv- 
ing the Fokker-Planck equation for a pure random walk; 
it leads to a term proportional to the Laplacian of P a in 
velocity space. The third term can be evaluated by insert- 
ing the expression (|178[) in Eq. (|185[) . Finally, we obtain 
the exact equation 



dPg 

dt 



dP a 
dr 



xP„^(r,v;r',t) = — 



j dr' ^2 mpVuQr - r'|) 

(186) 



d 


/* 


K 


0P a 


dv 



iaPaV 



where the statistical correlations are encapsulated in the 
two-body probability distribution 

P a ,p(r, v; r', t) = (5(r - r a (t))8(v - v a (t))5(r' - r„(t))). 

(187) 
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If we implement a mean field approximation 

P a ,p (r , v; r',t) = P Q (r , v, t)P (r' , t) , 
the preceding equation can be rewritten 



Substituting Eq. (TT!)g]) in Eq. pM]) , we obtain 



(188) 



dp e 



N 



dt 



^mV ■ (V$ ex {r a (t),t)8(r - r a {t))) 



5P Q 5P tt g>P q 

~m +v 



dr dv 



V J dr'Y,™ u(\r-v'\)Pe(r',t) =~V ■ (v$ ex (r,t)^m6(r-r a (t))j ,(197) 

so that, finally, 

dpex 



0_ 
dv 



D 



, dP* 
a dv 



£aP*V 



(189) 



If we extend the sum over all (3 and introduce the spatial 
density (163|) . we obtain 



dt 



(198) 



<9f 



<9r 9v 



(190) 



This exact equation, where p ex is expressed in terms of 
Dirac functions and <P ex is given by Eq. (|194[) . contains 
exactly the same information as the deterministic equa- 
tions (155)1 . 

For the inertial model, at T — 0, the equations of mo- 
tion are 



This can be rewritten in the form of a mean field Kramers 
equation 



8P a , y .^_ w dP c 



dt 



Or 



d 



dv dv 



D 



, dP» 
a dv 



C«PqV , 

(191) 
(192) 



dt 



dr a 

dt 



-£v Q - V$ ex (v a ,t). 



(199) 
(200) 



where 

<P(r,i) = y /)(r',t)u(|r- r'|)dr', 
is the potential. 

B The case T = 



Here, we assume that all the particles have the same mass 
m. At T = (no noise), the equations of motion (jTUJ) in 
the overdamped limit become 



Introducing the exact distribution function 



JV 



/«x(r,v,t) = ^m£(r-r Q (i))5(v-v a (t)), (201) 

a=l 

and taking the time derivative of Eq. (|201[) , we get 

N 



dfea 

dt 



mV • (r a <y(r - r a (t))5(v - v a (t))) 



Q = l 



JV 



]T m— • (v a 5(r - r a (t))6(v - v a (t))) . 



(202) 



(193) 



where <? ex (r,t) is the exact gravitational potential that is 
solution of the Poisson equation 



Substituting Eqs. $MB~ (2001) in Eq. (|2T)2"j) and following a 
procedure similar to that developed previously, we obtain 



Of, 



A<P ex = S d Gp ex , 
with the exact density field 

TV 

Pex(r,t) = Y mS ( r ~ r a(t))- 
a=l 

Taking the time derivative of Eq. (|195p . we get 

d N 

^ = -X)mV-M(r-r (t))). 



(194) 



(195) 



dt 



ex 

V • — V$ex 



dv 



dv 



f|;(/«v). (203) 



This exact equation, where f ex is expressed in terms of 
Dirac functions, contains exactly the same information as 
the deterministic equations (|199p - (|200|) . For £ = 0, it re- 
duces to the Klimontovich equation of plasma physics. 



C The case N = 2 in one dimension 

In this Appendix, we study the statistical mechanics of 
(196) t w0 particles in gravitational interaction in d = 1 where 
explicit results can be obtained. 
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C.l The microcanonical ensemble 

In d = 1, the Hamiltonian of two particles in gravitational 
interaction can be written 



10 



H = ~MV 2 + i/V + Gm im2 \r\, 



(204) 



where we have used the coordinates (R, V) of the center 
of mass and the coordinates (r, v) of the reduced particle. 
We note that the energy takes only positive values. Intro- 
ducing xi = {M/2) 1 / 2 V and x 2 — (fi/2) 1 / 2 v, the surface 
of phase space with energy less than E can be rewritten 



r(E) 



8R 
/Mjl 



dr 



x ll^ V E— Gra\iriir 



d 2 x. (205) 



The last integral represents the surface of a disc with ra- 
dius \JE — Gm\m 2 r. Thus, we get 



r(E) 



8ttR 



(E — Gm 1 m 2 r)dr. (206) 



The range of integration r m is such that E — Gm\m, 2 r > 
and r < R. Therefore, r m — E j Gm\m 2 if E < Gm\m 2 R 
and r m = R if E > Gm,\m 2 R. The density of states in the 
microcanonical ensemble is g(E) — dr/dE. Introducing 



A = 



8tt 

[M\L 



E 



we get 



Gm\m 2 R 



g(e) = ARr ri 



t = 



Gm\m2R 



(207) 



(208) 



with r m = eR if < e < 1 and r m = R if e > 1. 

The caloric curve t = t(e) and the equation of state in 
the microcanonical ensemble can be obtained from Eqs. 
(I133p - ()135|) and the exact expression <|208|) of the density 
of states. For < e < 1, we obtain 



5 (e) = AR 2 e, 



k B T = E -► t = e, 



PV = 2k B T. 



For e > 1, we obtain 



g(e) = AR 2 7 



k B T = +oo — > t = +oo, 



PV = 2k B T = +oo. 



(209) 



(210) 



(211) 



(212) 



(213) 



(214) 



The caloric curve in the microcanonical ensemble is rep- 
resented in Fig. [U 



6 - 



2 - 




Fig. 4. Caloric curve in microcanonical (MCE) and canoni- 
cal (CE) ensembles for the binary star model N = 2 in one 
dimension. The caloric curves are similar in the two ensem- 
bles (without negative specific heat region) so that there is no 
phase transition contrary to the situation in d = 3 [415] . There 
exists equilibrium states for all accessible energies E > in 
MCE and for all temperatures in CE. It can be of interest to 
compare these curves with the caloric curve obtained within 
the mean field approximation valid for N 3> 1 which presents 
a similar behaviour [12] . 

C.2 The canonical ensemble 

For N = 2 particles in gravitational interaction, the par- 
tition function in the canonical ensemble can be written 



Z 



e -W— e -P^ ■ e -f }mim ' G WdRdVdrdv. 



(215) 



Integrating over the velocities and over the position of the 
center of mastH, we get 



Z = 



8ttR 



-f3Gm 1 m 2 r i 



Therefore, the partition function is given by 
8ttR 



1 



-(3Grn\7Ti2B- 



(216) 



(217) 



/3 2 G(mim 2 ) 3 / 2 
The pressure can be computed from Eq. ([52")) leading 



to 



PV = k R T 



1 + 



Gm\ra2(iR 



(218) 



For T -> +oo, we have PV = 2k B T and for T -> 0, we 
have PV = k B T (this is the pressure created by an effec- 
tive single particle resulting from the collapse of the two 



In fact, the domain of integration of r depends on R. There- 
fore, our treatment (see also Sec.[5)l is approximate and will be 
improved in a future work. 
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Fig. 5. Distribution of energies for the binary star model 
N = 2 in one dimension in the canonical ensemble. We have 
plotted this distribution for different values of the temperature. 

particles). The average energy obtained from Eqs. (|15ip 
and (|2T7|> is 



(e) = 2t- 




(219) 



It behaves like (e) ~ 2t for t — > and like (e) ~ t for 
t — ► +oo. The caloric curve in the canonical ensemble is 
represented in Fig. |4] and it is compared with the micro- 
canonical caloric curve. Finally, the energy distribution at 
temperature T in the canonical ensemble obtained from 
Eqs. Ipgjl . ([21~7f and pOS]) is given by 

p{e) = *{i-e-v*y (0 - 6 - 1} ' (220) 

P ^ = t»(l - e-V«) ^ ^ 

The distribution of energies is represented in Fig. [5] for 
different temperatures. 



References 

1. Dynamics and Thermodynamics of Systems with Long 
Range Interactions, edited by T. Dauxois, S. Ruffo, E. Ari- 
mondo and M. Wilkens, Lect. Not. in Phys. 602, Springer 
(2002). 

2. P.H. Chavanis, Physica A 361, 55 (2006). 

3. P.H. Chavanis, Physica A 361, 81 (2006). 

4. T. Padmanabhan, Phys. Rep. 188, 285 (1990). 

5. P.H. Chavanis, Int. J. Mod. Phys. B 20, 3113 (2006). 

6. D. Lynden-Bell and R. Wood, Mon. Not. R. Astron. Soc. 
138, 495 (1968). 

7. P.H. Chavanis, A&A 381, 340 (2002). 

8. J. Binney and S. Tremaine, Galactic Dynamics (Princeton 
Series in Astrophysics, 1987). 



9. D. Heggie and P. Hut, The Gravitational Million-Body 
Problem: A Multidisciplinary Approach to Star Cluster Dy- 
namics (Cambridge University Press, 2003). 

10. W.C. Saslaw, Gravitational Physics of Stellar and Galactic 
Systems (Cambridge University Press, 1985). 

11. P.H. Chavanis, C. Rosier and C. Sire, Phys. Rev. E 66, 
036105 (2002). 

12. C. Sire and P.H. Chavanis, Phys. Rev. E 66, 046133 (2002). 

13. P.H. Chavanis and C. Sire, Phys. Rev. E 69, 016116 
(2004). 

14. C. Sire and P.H. Chavanis, Phys. Rev. E 69, 066109 (2004). 

15. P.H. Chavanis and C. Sire, Phys. Rev. E 70, 026115 (2004). 

16. C. Sire and P.H. Chavanis, Banach Center Publ. 66, 287 
(2004). 

17. P.H. Chavanis, M. Ribot, C. Rosier, C. Sire, Banach Center 
Publ. 66, 103 (2004). 

18. J. Sopik, C. Sire and P.H. Chavanis, Phys. Rev. E 72, 
026105 (2005). 

19. P.H. Chavanis and C. Sire, Phys. Rev. E 73, 066103 
(2006). 

20. P.H. Chavanis and C. Sire, Phys. Rev. E 73, 066104 
(2006). 

21. P.H. Chavanis, C. R. Physique, 7, 331 (2006). 

22. P.H. Chavanis, A&A 356, 1089 (2000). 

23. P.H. Chavanis, Eur. Phys. J. B 52, 433 (2006). 

24. P.H. Chavanis, Statistical Mechanics of Two-dimensional 
Vortices and Stellar Systems in Dynamics and Thermody- 
namics of Systems with Long Range Interactions, edited by 
T. Dauxois, S. Ruffo, E. Arimondo and M. Wilkens, Lect. 
Not. in Phys. 602, Springer (2002) |cond-mat/0212223| . 

25. P.H. Chavanis, Eur. Phys. J. B 54, 525 (2006). 

26. J. Sopik, C. Sire and P.H. Chavanis, Phys. Rev. E 74, 
011112 (2006). 

27. W. Gotze, Liquids, Freezing amd the Glass Transition (Les 
Houchcs) (North-Holland, Amsterdam, 1989). 

28. W.B. Russel, D.A. Saville, and W.R. Schowalter, Col- 
loidal Dispersions (Cambridge University Press, Cam- 
bridge, England, 1989). 

29. D.S. Dean and A. Lefevre, Phys. Rev. E 69, 061111 (2004). 

30. A. Einstein, Investigations on the Theory of Brownian 
Movement, edited by R. Fiirth and translated by A.D. 
Cowper (Dover, New York, 1956). 

31. A.M. Salzberg and S. Prager, J. Chem. Phys. Lett. 38, 
2587 (1964). 

32. A.M. Salzberg, J. Math. Phys. 6, 158 (1965). 

33. J. Katz and D. Lynden-Bell, Mon. Not. R. Astron. Soc. 
184, 709 (1978). 

34. M. Kiessling, Comm. Pure Appl. Math. 46, 27 (1993). 

35. E. Caglioti, P.L. Lions, C. Marchioro and M. Pulvirenti, 
Commun. Math. Phys. 143, 501 (1992). 

36. A. Alastuey, Ann. Phys. Fr. 11, 653 (1986). 

37. J.D. Murray, Mathematical Biology (Springer, Berlin, 
1991). 

38. E.F. Keller and L.A. Segel, J. Theor. Biol. 30, 225 (1971). 

39. T.J. Newman and R. Grima, Phys. Rev. E 70, 051916 
(2004). 

40. P.H. Chavanis, C. Sire, Physica A (in press) 

41. W. Jager, S. Luckhaus, Trans. Amer. Math. Soc. 329, 819 
(1992). 

42. P.H. Chavanis, Physica A (in press) 

43. G. Kirchhoff, in Lectures in Mathematical Physics, Me- 
chanics (Teubner, Leipzig, 1877). 

44. L. Onsager, Nuovo Cimento, Suppl. 6, 279 (1949). 



